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We investigate systematically sample-to sample fluctu- 
ations of the probability r of no return into a given en- 
trance channel for wave scattering from disordered systems. 
For zero-dimensional ("quantum chaotic") and quasi one- 
dimensional systems with broken time-reversal invariance we 
derive explicit formulas for the distribution of r, and inves- 
tigate particular cases. Finally, relating r to violation of S- 
matrix unitarity induced by internal dissipation, we use the 
same quantity to identify the Anderson delocalisation tran- 
sition as the phenomenon of spontaneous breakdown of S- 
matrix unitarity. 



Various aspects of chaotic wave scattering in the pres- 
ence of absorption or internal losses attracted a consid- 
erable attention in recent years [1]- [10]. In the gen- 
eral case a convenient framework for extracting univer- 
sal properties of the corresponding S-matrix is provided 
by the method of effective non-Hermitian Hamiltonian 
H = H — ITT X^ali W^aiSiM^a , m tcmis of which the energy- 
dependent element Sab of the AI x M scattering matrix 
S is expressed as: 

Sab=Sab-2t7rW^j^Wb, (1) 

see [11], [12] and references therein. Here H stands for 
a self- adjoint Hamiltonian describing the closed counter- 
part of the disordered or chaotic system under consid- 
eration, £ stands for the energy of incoming waves and 
the (energy-independent) vectors Wa, a = 1,2, ...M con- 
tain matrix elements coupling the internal motion to one 
of open M channels. As is easy to verify such a con- 
struction ensures the unitarity of the scattering matrix: 
S^S = 1m, provided the energy £ takes only real values. 
When one allows for energy parameter to have nonzero 
imaginary part e = lm£ > 0, the 5— matrix unitarity is 
immediately lost: S''' S < 1m. Physically the parameter 
e stands for a uniform damping inside the system, and 
is responsible for the fact of losses of the outgoing flux 
of the particles as compared to the incoming flux. The 
balance between the two fluxes is precisely the physical 
mechanism standing behind the S— matrix unitarity. 

In the present paper we concentrate on the "probabil- 
ity of no return" (PNR), which is defined as the quantum- 
mechanical probability for a particle entering the system 
via a given channel a never exit through the same chan- 
nel. This quantity is well-defined for a given realization 



of disorder and will show sample-to-sample fluctuations 
whose statistics we are going to study. In the case of no 
internal dissipation PNR is the same as the probability 
to exit via any of the remaining channels, known as the 
transmission probabihty J2b^a l^ab]"^ = 1 — [S'aal^. For a 
system with absorption the last equality is violated, and 
we keep the notation = 1 — [S'aaj^ for the PNR related 
to the reflection probability Ra in the same channel as 
Tq = 1 — Ra. In particular, if only a single open channel 
is attached to our disordered system and the boundaries 
are purely reflecting, then neglecting dissipation trivially 
results in T = 0. The nontrivial statistics of r then arises 
solely due to an absorption, and for small absorption the 
PNR value r ~ 2eTw [2,3,7], where tw is the so-called 
Wigner delay time intensively studied in recent years, see 
[3,4], [11], [13]- [15] and references therein. 

It is convenient to write explicitly the normalisation 
of the channel vectors as W^Wa = 7a/7r, and assume 
that different channel vectors are orthogonal: W^Wb = 
for a b. In fact, one should remember that the ef- 
fective strength of every open channel is more appropri- 
ately characterised by the so-called "transmission coeffi- 
cients" [12] (also known as the "sticking probabilities"): 
Ta = 1 — \ {Saa)\^, related to bare couplings -fa by 
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where v stands for the mean spectral density. Here and 
afterwards the angular brackets stand for the disorder- 
averaged value of the quantities. Two limiting cases 
Ta = 1 and Ta = correspond to situations of perfectly 
coupled and decoupled (closed) channel a, respectively. 

The starting point of our analysis is based on the 
following convenient representation for the diagonal el- 
ements of the scattering matrix: 
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where Ha = H-iir ^^^^ Wb(E)W^. In this way we reduce 
our problem to investigating statistics of the diagonal 
entries of the resolvent Ga of the "reduced-rank" non- 
Hermitian operator Ha independent of the vector Wa ■ In 
particular, for the single-channel case M — 1 the operator 
Ha will not contain the channel vector W at all, and will 
be therefore a self-adjoint one: H — H . 

The statistics of the diagonal entries of the resolvent 
of a random self-adjoint Hamiltonian H describing the 
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motion of a quantum particle in a static disorder was dis- 
cussed in much detail by Mirlin and Fyodorov [16] in the 
framework of the supermatrix nonlinear a— model [17]. 
In particular, for the case of systems with broken time 
reversal invariance they were able to find a very compact 
representation for the joint probability density V(u^ v) of 
the real u — KeGa and imaginary v = ImGa parts of 
the quantity Gq, assuming normalisation W\Wa = I/tt. 
Physically the variable v is the most important and be- 
ing the local density of states enjoyed thorough investi- 
gations, see [18], [19] and references therein. 

In fact, it is quite straightforward to incorporate non- 
Hermitian part of the Hamiltonian 7ia into the method, 
as was already partly done in [20] where statistics of 
ImGa was addressed as describing fluctuations of the 
photodissociation cross-section. 

According to [16] the function V{u,v) is given, for the 
centre of spectrum Ref = 0, by V{u,v) — j^P{x), 
where x stands for the combination x = {u'^ +v'^ + l)/2w, 
and the function P{x) is given by 



P{x) = LF{x), F(; 
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Hx,X), (4) 



where we introduced the (Legendre) operator L = 
^(.T^ — 1)^. The particular form of the function T{x, A) 
depends crucially on the effective spatial dimension of 
the underlying disordered system, and is, for example, 
quite different for zero-dimensional systems ("quantum 
chaos") and for "diffusive" extended quasi-one dimen- 
sional, or higher dimensional systems. We will give ex- 
plicit analysis of several physical possibilities later on in 
the paper. 

Having at our disposal the expression for P{u, v) we 
can relate the PNR distribution P{Ta) to the function 
P{x). After a set of algebraic transformations we find 
the following attractively simple formula: 
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x{Ta,e) = 1 — - 1 +4cOS^-5^^ '- 

\Ta / \Jo / TaJ-a 

Now we proceed with a separate analysis of a few physical 
situations possible in disordered systems. In all cases we 
assume time reversal symmetry to be broken. 

I. "Zero-dimensional" quantum chaotic system. 

We assume that the disordered region is coupled to M 
scattering channels characterized by effective coupling 
constants gi,...,gM, see eq.(2), with gi corresponding 
to the chosen entrance channel. The strength of uni- 
form damping will be characterized by the parameter 
7] = 27re/A, where A stands for the mean level spacing 
generated by the Hermitian Hamiltonian H. According 
to the standard argumentation, H can be effectively re- 
placed by a large N x N random Hermitian matrix taken 



from the Gaussian Unitary Ensemble, see e.g. [11,12]. 
Then in the limit of large enough N ^ M the function 
J^{x, A) depends on the remaining M — 1 coupling con- 
stants, as well as on the effective damping r] as [20]: 
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9a + X 
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The function F{x) in Eq.(4) can be found in a closed form 

for any M as one gets, in fact, a simple recursion relating 
Fm{x) to Fm-i{x). Here we restrict ourself mainly to the 
cases of one and two open channels M = 1,2: 



F,{x) = / — e-"-, F2(x) = - 2 

Jx—l 



s ''^ sinh 
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The distribution V{t) for M = 1 is then equal to (cf. [4]) 
2 

V{t) = -^e'""^ {loivB) [sinhr]ir]A- 1) -Ft? cosh??] 
- riBsinhrjhirjB)} , A = - l) (^^ - lY (8) 



B = 4 



V(l-r)(l-ri) 



where Io{z), Ii{z) stand for the modified Bessel functions 

of the respective order. For particular case of perfectly 
coupled channel Ti = 1 Eq.(8) reduces to the formula 

Vi{t) = ^e-^^'^ [t (1 + 2?7 - e^") + 2?? (e^" - l)] (9) 

derived earlier [2] with a very different method. 

The function V{t) for M = 2 can be obtained straight- 
forwardly, but the general formula is too long, and we 
restrict our discussion by a few particular cases. First 
of all, when dissipation is absent (77 = 0) we recover the 
exact distribution of the transmission probability found 
earlier in [21,22] by rather different methods. Next case 
to be considered is that of a lossy system coupled to two 
perfectly open channels 51 = 52 = 1: 



V2{t)=Vx{t) + 



1 - e-^''/^ 

■n 
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where Pi{t) is given in Eq.(9). In fact, it is not difficult 
to find a similar recursive formula relating P{t) for M 
perfectly coupled channels to the same function for Af — 1 
perfect channels. We do not give that formula, apart 
from the simplest case of no dissipation: 



Vm{t)^Vm-i{t)+[{M-1} 



.M-2 



(M - 2)t 



M-31 



M > 

M-3 



which immediately yields Vm{t) — {M — \)t^ ^ 
1. This formula (as well as its counterpart V{t) oc 
for preserved time reversal invariance), in fact, follows 
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from the known distribution of 1 x 1 subunitary block of 
random imitary scattering matrices, see [23]. 

II. Quasi- ID systems. Consider a single channel at- 
tached to one edge of a piece of quasi one-dimensional dis- 
ordered metal of length L, with the opposite edge being 
in contact with perfectly conducting lead of very many 
channels. When the internal dissipation is absent, the 
function F{x) was found by Mirlin [24]: 



F{x) = In 



X ■ 

dkk 
1 + 



(11) 



tanh 



where the dimensionless parameter t = L/S, is the sam- 
ple length L measured in units of the localisation lengths 
^. The (real) functions Pi,{x), v = —l/2 + ik/2 are 
known as conical functions, and represent a special case 
of Legendre functions. As such they satisfy: IjPv{x) = 
v{v + \)Piy{x). This observation immediately yields the 
following expression for the PNR distribution: 



(12) 



I{t;x)= dkktanh(^^^P_^^iM{x)e-*^'''+^'>/\ (13) 

where we assumed, for simplicity, that the selected single 

channel is perfectly coupled to the scattering medium. 
Surprisingly, this distribution is practically the same as 
the distribution of the reflection coefficient from a piece 
of strictly one-dimensional medium obtained long ago in 
the framework of the Berezinskii technique [26]. In par- 
ticular, for any value of the parameter t the distribution 
displays a log-normal far tail corresponding to very small 
PNR values r — > 0. To find it for t <C 1 one needs an 
asymptotic of the conical functions for large arguments 
a; 1 which we borrow from Eq.(50) of the paper [24]. 
A calculation very similar to that presented in [24] yields: 



.X 1 V— Inr 
V{t < 1) ~ —f^- exp 
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t < |lnr| 
(14) 



This log-normal tail is related to the presence of the 
anomalously localised states [25]. In the opposite case 
of very long samples t 1 the PNR values are expo- 
nentially small due to localisation phenomenon and the 
distribution is purely log-normal: 
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1 1 
2^7^ 



exp 



-l(. + ln.)^ 



t ■ 



|lnr| > 1 
(15) 



Let us turn our attention now to the case of a quasi- ID 
disordered sample with a nonvanishing internal dissipa- 
tion e > 0, assuming second edge of the sample to be 



impenetrable for waves. The scaling physical parameter 
controlling the role of dissipation is then given by [18] 
S = npe^, with p standing for the mean spectral density. 
This is just the dissipation e measured in units of the 
mean level spacing for a sample whose length is ^. The 
most interesting regime is that of small 5 <C 1. In that 
limit the function J^{x, A) turned out to be independent 
of A, whereas the x— dependence persists in a form of 
the scaling combination y = 2Sx, i.e. J^(x,\) F{2Sx). 
This implies that the relevant values of parameters are 
2/r~a;'~(5~^ » 1. The latter condition immediately 
results in the formula F{x) A5F{y) /y, and also con- 
verts the Legendre operator L to Ly = ^j/^^. Let us 
note that the emerging PNR distribution yields, in fact, 
the distribution of the Wigner delay time via the relation 
Tw = 27r/9^/t/. 

The expression for the function F{y) is known explic- 
itly [27]: F{y) = F^{y) + Ft{y), where 



Ft{y) 



dkk Trfc 



(2y^) 



-t(fe^ + l)/4 



(16) 



and Faciy) = 2y^Ki {2y^^, with K^{u) standing for the 

Macdonald function. 

For the case of very short (t <C 1) sample the func- 
tion F{y) is known to be approximated by exp — (fy) 
[27]. A simple calculation then yields the distribution 
V{t) = [{U5f/T^] exp{-4M/T}. Realising that 2t5=ri 
we see that the distribution coincides with the weak ab- 
sorption limit of Eq.(9). As expected, the same distribu- 
tion follows from that of the Wigner delay time [11]. 

In the opposite limit of very long samples t — > oo only 



first term survives, and by noticing that J-iy 



Foo{y)/y 



Poo{y), we find the corresponding PNR distribution: 



V{t) = 



16 /5 



5/2 



Ki (4VV^) . (17) 



Although the typical value of r is of the order of 6, the 

moments (r™) do not exist for m > 1 because of the 
power law tail V{t ^ S) oi t~^. A similar tail was found 
in the distribution of the total reflection coefficient from 
multichannel long disordered ID sample in [9] , and is also 
typical for the Wigner delay time distribution in purely 
ID system [15]. Negative moments of t are equal to 
(r-*^) = (45)-'=fc!(A: + 1)!. Note that they differ from the 
corresponding moments in purely ID system [15] by extra 
factorial factor {k + 1)!, reminiscent of similar relations 
between other quantities in ID and quasi ID [18]. 

Finally, in the case of strong absorption 5 » 1 in a long 
wire t ^ oo the function J^{x, A) = cxp{—\/5{x — \)} [17] 
and the resulting distribution 'P(r) coincides with that 
given by Eq.(9), with r] replaced by y/S. 

III. Behaviour at the Anderson transition. 

Lets us shortly discuss a possible qualitative behaviour 
of the PNR T in a scattering system formed by a single 
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perfect channel attached to a d— dimensional disordered 
sample at the vicinity of the point of the Anderson de- 
localisation transition ac (the mobility edge). Here we 
denote by a an effective parameter which controls the 
transition in the infinite sample, with states being lo- 
calised (extended) for a > Uc (respectively, a < ac)- 

Our arguments are based on a picture of the transition 
as described in terms of a functional order parameter de- 
veloped in detail in [16], sec also earlier results in [17] and 
[28]. For a sample of finite size L the PNR is a function 
of three parameters: e,L,a. According to the suggested 
scenario, the behaviour of the function J^{x, A) in the in- 
sulating phase in the weak absorption limit 5 cx e^"^ — > 
is expected to be reminiscent of that described above for 
the one-dimensional case, i.e. .F(a;, A) F{2Sx), and the 
function F(y) decays to zero for j/ 3> 1. Then it is natural 
to expect that all negative PNR moments in the infinite 
volume limit L ^ oo are to scale as (t^*^) ~ ^-k^-dk^ 
where ^ is the localisation length diverging in the vicinity 
of the mobility edge. 

In contrast, in the delocalized phase the function 
J^(.T, A) is expected to remain a non-trivial function of 
both X and A even when e ^ 0, provided the latter limit 
is taken after the infinite volume limit L oo. This 
should immediately result in a finite-width distribution 
P{t) of the PNR. From this point of view the Ander- 
son transition acquires a natural interpretation as the 
phenomenon of spontaneous breakdown of S-matrix uni- 
tarity. As long as a ^ ac, the widths of the distribution 
and properly defined (negative) PNR moments should 
vanish, with some set of critical exponents. 

If, however, we take limit e ^ first, then for a < ac 
PNR in large but finite sample should scale with the 
system size L as (r"*^) ~ C{a)e~^ L~'^^ , where C{a) 
is expected to diverge when a — » ag. In some sense 
the behaviour of the negative moments of the Wigner 
delay time defined as tw = ^va\^^QT{e,L,a)/2e is rem- 
iniscent of that for the inverse participation ratio [25]. 
This analogy suggests a possibility for anomalous scaling 
{t^) ^ L~'^^'' with ffe ^ A; at the mobility edge a = ac, 
which would then reflect the underlying multifractality 
of the wavefunctions. 

It will be very interesting to perform a detailed numer- 
ical analysis of PNR and Wigner delay times for realistic 
and well-controlled models of scattering from disordered 
systems, e.g. quantum graphs [29] or models used in 
[13] and to verify the suggested picture qualitatively and 
quantitatively in various regimes. The statistics of PNR 
should be also of experimental accessibility in microwave 
resonators type of experiments, see e.g [6,7] and refer- 
ences therein. 
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